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ABSTRACT 


The objectives of the project described in this 
thesis were to : 1) provide the documentation that is 
needed for a Naval Postgraduate School student to use 
the general purpose finite element computer progran 
called SAP IV, and 2) to make available, and prepare 
the users manual for, a pre-and postprocessor progran 
called SUBROUTINE PSAP. This subroutine, which was 
developed at the NASA Langley Research Center, has 
been modified to specifically plot the finite element 
modei gecmetry for SAP IV models and to postprocess 
displacement data for those models on the NPS Calcomp 
Model 765 Plotter. The input and output for SAP IV and 

SUBROUTINE PSAP are discussed in detail. The codes 
have been used successfully in AE 4102, Plight Vehicle 


Structural Analysis II. 
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A. STRUCTURAL ANALYSIS SOFTWARE PACKAGES 


The past decade has seen great strides in the field of 
computer software packages that have been developed for use 
in structural analysis. A great deal of time and effort has 
peen expended in the development and distribution of these 
packages. Today there are structural analysis computer 
programs for almost every conceivable structure an engineer 
could desire to analyze. The types of analyses performed by 
a given package vary widely and may include features for 
linear or nenlinear materials, static analysis, dynamic 
analysis, buckling analysis, and nonlinear dynamic analysis, 
to name only a few. These programs can be grouped into two 


major categories - special purpose or general purpose. 


In order to expose the students at the Naval 
Postgraduate School to the use of a general purpose 
structurai analysis program, as well as to provide the 
capability of using such a program in research work, the 
finite element structural analysis program SAP IV (Reference 
1) was acquired and made operational at the Naval 
Postgraduate School by Professor Gilles Cantin of the 
Mechanical Engineering Department. SAP IV can perforn 


linear static and dynamic analyses on one-, two-, and 


three-dimensional structures. 


B. PRE-AND POSTPROCESSORS FOR STRUCTURAL ANALYSIS PROGRAMS 


After the development of a large number of structural 
analysis programs, users began to recognize that a disparity 
existed between efficient general purpose structural 
analysis programs and optimum utilization of these programs. 
Many of today's software structural analysis packages, such 
as SAP IV, require the user to prepare and reduce tremendous 
amounts of data. The need existed for some aids in 
processing and reducing these large quantities of data. 
Consequently, there have been many pre- and postprocessors 
developed for a specific use as well as for a general use 
basis during the past several years. The value of a given 
processing fpackage lies in its ability to aid the user in 
preparing his model geometry, in his data checks, and in 
processing the output in an easily understood fashion. One 
of the most effective means utilized in preparing or 
reducing data is through the use of visual displays, whether 
they be designed primarily for graphic presentations, such 
as the Strcmberg-Carlson, or for paper plots, such as 
Calcomp. This ability to visually display input and output 
data is a highly valuable tool for the structural analyst. 


Cs THESIS MOTIVATION 


The desire to hawe a pre- and postprocessor that could 
be used in conjunction with the general purpose structural 
analysis program SAP IV prompted the acquisition and 
implementaticn of a genegal use plotting package by this 
author. After researching the possible options, an existing 
program was obtained from Anamet Laboratories, San Carlos, 


California. The program was originally developed at the 
Langley kesearch Center, Hampton, Virginia by Gary L. Giles 
ror use with modern digital computers. The program, details 
of which can be found in Reference 2, generates oblique 
orthographic projections of three-dimensional finite element 
nodels and is distinguished by its provisions for 
generality, ease of use, different display options, and 
computational speed. The computer code was written for use 
on CDC 6000 series machines and had to be modified somewhat 
for use on the NPS 360/67. The modified version of the 
program is now available for use with the NPS Calcomp 
Plotter Mcdel 765. 


The primary purpose of this thesis is to provide the 
hecessary documentation in order that students enrolled at 
the Naval Pcstgraduate School, and specifically in the 
courses AE 4101, 4102 (Flight Vehicle Structural Analysis I 
and II respectively ), can, with a minimum of difficulty, 
effectively utilize SAP IV and its now-functional pre-and 
postprocessor PSAP. The remainder of this thesis is broken 
down into two major subdivisions, 

1) IL. GUIDE TO THE USE OF SAP IV AT THE NAVA 
POSTGRADUATE SCHOOL 


== SS eS ee ee 


2) III. GUIDE T2 THZ USE QF SUBROUNTINE PSAP, A 
PRE-AND FOSTPROCESSOR FOR SAP IV 
Appendices A and B of this work provide detailed examples 


for the input preparation and output reduction of data using 
both SAP IV and SUBROUTINE PSAP. 
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aA. DESCRIPTION 


SAP IV is a general purpose structural analysis digital 


computer program that can provide a finite element solution 


for both the static and dynamic analysis of linear 
st 


tructural systems. A detailed user's manual is contained 
in Reference 1. The program has the capacity to analyze very 
large three-dimensional systems, as well as small systems, 
with no loss in efficiency. SAP IV, which is coded in 
FORTRAN IV , is a very flexible program and can be 
considered a very efficient aid to the analyst. The purpose 
of this section is to provide the necessary additional 
documentation, above that provided in Reference 1, for a 
Student at the Naval Postgraduate School to make use of the 


progran. 


The metaods of analysis and the construction of the 
program are not included in this section, but can be found 
in Reference 1. The program contains nine finite elements of 
the following types: 

(a) three-dimensional truss element, 

(b) three-dimensional beam element, 

(c) plane stress and plane strain elenent, 

(d) two-dimensional solid element, 

(¢) three- dimensional solid element, 

(f) variable-number-nodes thick shell and 
three-dimensional element, 

(g) thin plate or thin shell element, 


(h) boundary 2lement, 

(1) pipe element (tangent and bend). 
There ar2 numerous options and combinations of static and 
dynamic analysis that are available to the user of the 
program. Reference 1 provides specific details of the many 
available user options. 


B. COMPUTER CARD DECK PREPARATION AT NPS FOR SAP IV 


Figure 1 outlines the overall computer card deck 
necessary to access and utilize SAP IV as it is currently 
operational at NPS. A complete detailed breakdown of the 


necessary IBM job control cards follows on the next pade. 


4/ { STANDARD GREEN JB CARD ) 

4/GQ EXEC PGM=SAP,REGION=260K 

//STEPLIB OD UNIT=2321sVOL=SER=CELOO2,DISP=SHRy 
// OSN=FO0559.SAPLM 

//GO-FTCLFOOL OD UNIT=SYSDA,SPACE=(CY¥L2(696)),4 
¢/ OC B=(RECFM=VS3,LRECL=1284,BLKSIZE=2572) 
//GO-«FTO2ZFOO0!L OD UNIT=SYSDA,SPACE=(CYL2(606)),% 
/f CC B=(RECFM=VSB,LRECLHAL284, 8BLKS IZE=2572) 
//GOeFTC3FOOL DD UNIT=SYSDAy,SPACE=(CYL,(6,6)), 
4/ CC B=(RECFM=VSB,yLRECL=1284, BLKSIZE=2572) 
//GO«FTC4FOOL OD UNIT=SYSDA,SPACE=(CYL1(646)), 
// OC B=(RECFM=VS8,LRECL=1284, BLKSIZE=2572) 
4/GO-eFTCSFOOL DO ODNAME=SYSIN 


//GCeFTC6FOOL DD SYSOUTH=A,SPACE=(CYL, (351)), 
// OC B=(RECFM=FRA,LRECL=L33,BLKSIZE=1330) 
//GGeFTC7TFOOL OD UNIT=SYSDA,SPACE=(CYL4(6,5)), 
// OC B= (REC FM=VS3,LRECL=1L284, 8LKSIZE=2572) 
//GOeFTC8FOOL DD UNIT=SYSDAy_SPACE=(CYL (595) )y 
4/ CC B=(RECFM=VS8,LRECL=AL234, BLKSIZE=2572) 
//GOeFTCIFOOL DD UNIT=SY S04, SPACE=(CYL2(646)), 
// OC E=(RECFM=VS8,LRECL=1284, BLKSIZE=2572) 
//GOeFTLOFOOL DD UNIT=SYSDA,ySPACEH=HICYL 1 (676)), 
// OC B=(RECFM=VS38,LRECL=H1284, BLKSIZE=2572) 
//GQeFTLIFOOL DD SYSOUT=B,SPACE=(TRKy (2022) )¢ 
// OC B=(RECFM=F8, LRECL=80, 8LKSI ZE=7200) 
J/G9eSYSIN DD * 


€ SAG LV === “DATA I} 


/ ( STANDARC APS EOF CARD ) 
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/* 


| SEITE Elia CT ae ES eS 
ELEMENT LORD 
MULTIPLIES 


CONCENTRATED LORD 
MASS DATA 


ELEMENT DATA CARDS 
SAP IV DATA 


NOORL FOIN 


T DATA 
C 


// STANDARD GREEN JOS 
CRRD 


Figure i - SAP IV DECK SET-UP 
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After the //GO.SYSIN DD * card, the deck of SAP IV data that 
is prepared according to tae instructions in the appendices 
of Reference 1 follows. Following the SAP IV data is the 
standard NPS end-of-file card (/*). 


Ge HELPFUL POINTS ON DATA DECK PREPARATION 


There are several possible areas where a user can err in 
his preparation of data for SAP IV. The following items are 
representative of a few common sources of error. 

1) Particular attention should be directed toward 
uSing the correct formats and correct card columns in key 
punching data on cards. (i.e., integer formats, rigat 
justiried) 

2) The program has several internal data generation 
features inherent in it. Simply because data were generated 
during the program execution does not necessarily mean they 
were generated correctly. Errors in the user- prepared input 
cards used in data generation can cause severe discrepancies 
to occur during program execution. The generated data 
should be carefully checked for accuracy. 

3) In order to terminate a given problem a number of 
blank cards are necessary at the end of the input data deck. 
Specific points of interest pertaining to the problem of 
termination can be found in Reference 1 under note (1) on 
page II.1 and in note (1) on the top of page V.2. 


SAP IV has a wide range of options where a clever user 
can fully exploit the full capabilities of the program. The 
detailed decription of these features is found in Reference 
1; however, the following list summarizes a few of the more 
useful options and features. 


uk 


aero esa CS RA aE EN ERE PNA NAH Poe ee 


i a lh Sata et 


1) Data check only mode of execution. 
2) Nodal point and element data generation 
capabilities. 
3) Five different types of analysis. 
a. Static analysis 
b. eigenvalue /vector solution 
c. forced dynamic response by mode superposition 
d. cvresponse spectrum analysis 
e. direct step py step integration 
4) Automatic punched computer card output of 
displacement results. This feature was added to the progran 


by this autncr to make possible graphic postprocessing. 


The foregoing discussion is not intended to be a 
complete diagnostic summary of SAP IV, but ratner an aid to 
the student who desires to get started using the program as 
it currently exists on the NPS IBM 360/67. The above 
discussion, along with a copy of Reference 1, should provide 
a jumping off place for a novice structural analyst. A 
complete example of the job control cards and input data 


deck for a static truss analysis can be found in Appendix A. 


D. ALTERATION OF SAP IV AT NPS 


Because SAP IV is a very flexible program, it is 
possible for a user to make modifications to the basic 
program. The computer code is complex, but it is not overly 
dizficult to modify parts of the program in order to satisfy 
a specific user need. A method that has proven successful 
tor this author, in the modification of a SAP IV. subroutine 
to provide punched output of displacement data, and to 
create a personal version of SAP IV, is outlined in this 
section. 
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The first step is to define the subroutine name that the 
user wishes to modify. A listing of the source program is 
available through Professor R. E. Ball, Department of 
Aeronautics, or Professor G. Cantin, Department of 
Mechanical Engineering. Having defined those portions of the 
pregram for which modification is desired, the next step is 
to obtain a punched copy of those desired routines by using 
the following rormat. 

// (STANDARD NPS JOB CARD) 
//SYSPRINT DD SYSOUT=A,SPACE=(TRK, (10,1) ) 
//STSUT1 DD DISP=SHK, UNIT=2321, VOL=SER=CELOO1, 
// DSN=F0099.SAPSR 
//SYSUT2 DD SYSOUT=B 
//SYSIN DD * 
PUNCH TYPORG=P0,MAXNAME=2 
MEMBER NAME=(SUBROUTINE NAME) 
/* 
NOTE: SUBROUTINE NAME is the name of the SAP IV subroutine 
reguired. 


With the desired routine decks now in hand, the 
necessary changes can be incorporated into the subroutine 
deck and an cbhject deck is then obtained as follows: 

7/ (STANDARD NPS JOB CARD) 

// EXEC FORTCD 

// FORT.SYSIN DD * 

(MODIFIED FORTRAN SOURCE DECK) 
/* 


The next and final step is to take the object deck 
obtained in the previous step and insert it into the proper 
position in the control cards that are illustrated on the 


following page. 
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There are many Options available in the creation of load 
module libraries and they are discussed in detail in 
Reference 4, sections II. and [III. The portion of the 
preceding example control card deck that would necessitate 
modification is the J/ULINK.SYSLMOD card, where 
FQ559.SAPLM (SAP) should be changed according the following 
format : 

General format: Lnnnn. anyname 
where Lis ;: 
S-student data set 
F-faculty data set 
onnn is : 
user number assigned 
anyname is : 
any umigue name assigned by the user (1-6 
characters in length with the first character alphabetic), 
and CzL002, should reflect an appropriate data cell with 
available space. 


The procedure outlined above wiil create a modified 
version or SAP IV on @ chosen data cell available to the 
user for a period of 90 days. The program can now be 
executed with appropriate modifications to the //STEPLIB 
card as discussed previously in section II. B. Any questions 
concerning the creation of load module libraries or their 
execution can be answered by any of the programing 


consultants on the first floor of Ingersoll Halil. 


A. SUBROUTINE DESCRIPTION 


SUBROUTINE PSAP is a modified version of the oblique 
orthographic projection program that is found in Appendix B 
of Reference 2. The program, Originally developed for use 
at the NASA Langley Research Center, Hampton, Virginia, 
required some changes so that it could be used in 
conjunction with SAP IV and the NPS Model 765 Calcomp 
Plotter. The original version of the plotting package 
allowed for various geometry and displacement data input 
options. However, the subroutine, as currently filed, is 
constructed strictly for use with the input and output of 
SAP IV. With some slight aodifications to the subroutine, it 
could be adapted to any number of different types of input 
geometry or displacement data decks. A method for 
accomplishing this is outlined in section III. Cc. 


The current capability for generating obligue 
orthographic projections of SAP IV finite element nodelis is 
limited tc the following types of elements: 

1) Type 1, Three-Dimensional Truss Elements 

2) Type 2, Three-Dimensional Beam Elements 

3) Type 3, Plane Stress Membrane Elements 

4) Type 4, Two-Dimensional Finite Elements 

5) Type 6, Plate and Shell Elements (Quadrilateral). 
The undeformed topology of the finite element model, useful 


in checking input data, as well as the displaced topology 


| 
t] 
/ 
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projection of the same model, can be obtained from PSAP. The 
subroutine contains many different options and permutations 
of those cptions, some of which are listed below. 
1) plots of models annotated with grid point numbers 
of element numbers 
2) plots of portions of models 
3) exploded plots of model sections (i.e., line 
elements coincident with the edges of triangular or 
quadrilateral elements may be difficult to single 


out; program provides a capability to separate 
elements so that their absence or presence is easily 
detectable) 

4) displacements superimposed on grid point 


coordinates of the undeformed structure 

5) displacements represented as vectors extending the 

undisplaced grid points. 

There are many more Combinations of options that are 
available and will be more specifically outlined in the 
remainder of this section. Appendix B details a complete 


input-output example of a pin-jointed truss. 


Be. SUBROUTINE PSAP USE 


1. General Set-up of Input Deck 


In general, the correct sequence of computer cards 
required to utilize SUBROUTINE PSAP is shown schematically 
in Figure 2 and consists of eight separate major groups as 
follows: 

1) a group of JCL cards and main program to allocate 
storage 
2) a single card containing title information 


3) Namelist OPTION containing values to determine if 


a BA ME ig inl 


proper storage allocation is available and specifying 
various program options 
4) a geometry deck (SAP IV data deck) containing grid 
points and connectivity of the model 
5) an optional single title card used to identify the 
deck of displacement data to be plotted 
6) a single card containing the value of the total 
number of SAP IV load cases and an optional scale 
factor 
7) tne deck of displacement data to be plotted 
(output of the execution of SAP IV-static analysis) 
8) Namelist PICT containing values to specify the 
type of plot desired and what information is to be 
included on the plots. 
By repeating parts of the pasic input data to the prograa, 
different plots of the same data can be generated. 


The sequence of cards on the following page depicts 
the necessary JCL to execute SUBROUTINE PSAP. 


// (STANCAROD GREEN J38 CARD) 
// EXEC FORTCLGP,REGIONeGO=180K 


//FCRTeSYSIN CD * 


( MAIN PROGRAM ) 


/ ( STANCARE NPS E9F CARD ) 
J/LINKeUSCD CC UNIT=H2321,VIL=ESER=CELONZ,DISP=SHR, 
// OSN=FO559ePSAPLM 
//LTNKeSYSIN CO *¥ 
INCLUDE USOD( PSAP) 
ENTRY MAIN 
/ ( STANCARC NPS EGF CARD ) 
//SGeFTLOFOOL OO UNIT=SYSOA, 
// SPACE=(CYL 9(391)), 
// SCB=(RECFM=VS,BLKSIZE=3520) 


//GOeSYSIN DC # 


( CATA --- STARTING WITH THE TITLE CARD --- S&e FIGURE 2) 


/ ( STANCARC NPS ESF CARDO ) 
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/* 
REARINDER OF INPUT DECK 
REPEATS PARTS DEFENDING 
ON THE VALUE OF Kove iN 
NAMELIST PICT 


(FROM SRP IY 
EXECUTION?) 


DATRG CONTROL 
CARD 


DISPLACENENT ORTR DECK 


OMIT IF ae 
| DCRSE=2 NTLC, SCALEF 


CASE ID CARD 


GEOMETRY INPUT DATA DECK ie 


| 
NAMELIST OPTION 
TITLE CARD 


(ror SRP IY) 
UCL CARDS AND MRIN 
CALLING PROERAM 
Pigure 2 - PSAP INPUT CARD SEQUENCE 


ee 


3. Main Progran 


The main program consists of three cards. It 
allocates the proper storage and calls SUBROUTINE PSAP. An 
example of a main program for a finite element model that 
contains 400 nodes is : 

DIMENSION ZZZ (2800), DISPD (5,3,400) 
CALL PSAP (222, 2800,DISPD,400) 
END 
The definition of the arguments used in calling SUBROUTINE 
PSAP (22Z, NZ, DISPD, NON) are: 
ZZZ-blankcommon array used to store nodal coordinates 
and displacements 
NZ-length of array ZZZ (NZ is determined by muitiplying 
the number of nodes in the model by seven, i.e., NZ=no. of 
nodes * 7) 
DISPD-a three dimensional working array used in 
subroutine DATA9 i 
NON-number of nodes in the nodel 
It is crucial to dimension ZZZ (7*NON) and DISPD (5,3,NON) 
correctly in the main program. Improper dimensioning of 
these arrayS can cause output errors that are not readily 
traceable. 


4. Title Card 


This card contains any desired alphanumeric 
information in columns 1 to 80. The title will appear at the 
beginning of the plots. 
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Namelists are a convenient means of inputting the 


j names of several parameters along with their corresponding 
values. For the NPS IBM 360/67, the format for using 
| namelists is as follows: 

1) card 1-&name-beginning in card column 2 where name 
is the name of the subject namelist 

2) succeeding cards- beginning in card column 2, the 
hames of the variables and their values separated by 
commas 

3) final card -&£ND-starting in card column 2. 


An example of a namelist format is shown below: 


ce rT a ah Rial I 


1 x-direction displacements to be input 


col 
| 
. SOPTION 
NNDEST = 400,NUDISP =1, 
PSIZE =8.0 
SEND 
Any or all farts of a defined namelist may be included, and 
each parameter may be specified in any order between the 
€éname card and the &END card. 
The description of the variable names in Namelist 
OPTION and their default values are contained in Reference 
2. They are given here to assist the user in data 
preparation. 
FORTRAN name - Default value - Description 
NNDEST - 1 The number of nodes(NON) as defined in the 
progran 
NUDISP - 1 0 x-direction displacements not to be input 
: 


NVDISP - 


NWDISP —- 


KaTA - 9 


NVALUS - 
IRESEQ - 


KPLOT - 


XSPACE - 


1 


1 


1 


0 
1 


5.0 


ates RD nt = 


0 y-direction displacements not to be input 
1 y-direction displacements to be input 

Q z-direction displacements not to be input 
1 z-direction displacements to be input 

( NOTE: when SAP IV displacement data is to 
be used, NUDISP=NVDISP=NWDISP=1; for no 
displacement data NUDISP=NVDISP=NWDISP=0) 
Specifies the subroutine and corresponding 
method of input for model geometry 

1 subroutine GEOM1, a user-supplied 
Subroutine 

2 subroutine GEOM2, a user-supplied 
Subroutine 

9 subroutine GEOM9, reads in grid points aad 
element data specifically from a SAP IV data 
deck 

Specifies the subroutine and corresponding 
method of input for displacement data 

1 subroutine DATA1, a user-supplied 
Subroutine 

< subroutine DATA2, a user-supplied 
subroutine 

9 subroutine DATAI, reads a punched output 
displacement deck from execution of SAP IV 
NOT INCORPORATED-ALLOW DEFAULT 

Grid point numbers are stored in the prograa 
from 1 to the total number of grid points 

0 no internal resequencing of grid points 
necessary; they are already in ascending 
order starting with 1 

1 resequence grid points from lowest grid 
point number to highest grid point number 
Specifies the type of output device to be 
used ( ALLOW DEFAULT ) 

Space between plots in the y-direction,in 


inches( see Pigure 3 for an explanation of 
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axis orientation ) 


PSIZE - 10.0 Paper size in x-direction, in inches( used in 
scaling of plots to insure this dimension is 
not exceeded) 

IDCASE - 0 0 no identirication card preceeds the deck 
of displacement values 
1 identification card preceeds the deck 


of displacement values 


This portion of the input deck contains the grid 
point locations and the element connectivity. The deck has 
one ofr the following forums, depending on the value of KGEOM 
in the Namelist OPTION. 

KGEOM = 9 
Calls subroutine GEOM9, which is constructed to read 
SAP IV geometry data. 


(a) When KGEOM is specified as 9, the complete input 
deck of SAP IV data, prepared according to Appendices I 
through I¥ of Reference 1, is placed after the 6&END card of 
Namelist CPTION. The portions of SAP iV data deck that 
involve load, mass, or dynamic analysis data are not part of 
the input geometry data to SUBROUTINE PSAP. Only the grid 
point locations and the element connectivity data are used 
to generate the orthographic projection of the model. 
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SECOND 
PLOT 


PEN 
POS!T!ION 
AT END OF * 

PLOT 


XSPACE 


tL Fen POSITION AT 
BEGINNING OF PLOT 


A +X 


Figure 3 - NPS CALCOMP PLOTTER ORIENTATION 


(b) It is possible to input only a portion of the 
finite element model for a data check. To do so, it is 
necessary tc modify SAP IV element control cards (described 
in detail in Reference 1, Appendix IV) for the desired 
element types to reflect the portion of the element 
connectivity cards that will be input. For example, if only 
connectivity for element numbers 15 through 50 of element 
type 1 (truss elements) of a SAP IV data deck are available 
for input, it is necessary to alter the element control card 
for the truss elements. [he field on the SAP IV control card 
that defines the total number of truss elements (card 
columns 6-10) would reflect the upper bound, in this 
case-50, and card columns 65-70 would reflect the lower 
bound, in this case~15. All nodal coordinates for the 
entire amcdel may be input, or only those that specifically 
define the portion of the finite element model to be 
plotted. In either case, the nodal coordinates that relate 
to element numbers 15-50 must be specified. Unknown results 
will occur when trying to plot an element whose node points 
are not specified. The above feature is valuable in a case 
where several different people are preparing different parts 
of a large data base for a SAP IV problem and desire to 
individually check their inputs grapnically for accuracy. 
KGeON = 1 

Calls subroutine GEOM1, which is prepared py the user 
to read geometry data from a program other than SAP IV. 
KGEOM = 2 

Calls subroutine GEOM2, which is prepared by the user 
to read geometry data from a program other than SAP IV. 


Use of KGEOM=1 or 2 requires modification of 
SUBROUTINE PSAP to fit the specific format of the user's 
input gecmetry data. A method for doing this will be 
discussed in paragarph C. of this section. 
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If ILCASE =0 is specified in Namelist OPTION, this 
card is omitted. The card, if present, contains any desired 
alphanumeric information in card columns 1-80 which will 
identify all displacements for a given case. For IDCASE = 1 
and SAP IV puncaed displacement data, a case identification 
card must appear before each Namelist PICT for every 
different ioad case that is plotted in addition to load case 
one. This is illustrated in Figure 4. A maximum number of 
five different load cases can be obtained from SAP IV. The 
case ID card information, if present, will appear before 
each jload case's DISPLACEMENT DATA TO BE PLOTTED section in 
the printed computer output from SUBROUTINE PSAP. This 


information does not appear on any Calcomp plots. 


| 


(NOTE: ONE 0 
SAP IY CONT 
LORD 


UCL CARDS 


Figure 4 - 


Se cLIST PICT 
(LORDCRSE NO. 2/KDDE=2) 
ATA DECK FROM 


roree 


AINS ALL CASE 1D CARD 
CASES) CLORDCRSE NO. 2) 


alate tilelen tent ani ymin as 


DISPLACEMENT DATA DECK 
(FROM SRP IV) 


CRSE id CARD 


GECMETRY DATA DECK (SAP IY) 


SAMPLE SUBROUTINE PSAP EXECUTION DECK WITH 
IDCASE=1 


This deck contains the displacement values of the 
nodal points. The deck has one of the following forms, 
depending on the value of KDATA specified in Namelist 
OPTION. 

KDATA = 9 
Calls subroutine DAIA9, which reads SAP IV displacement 
data. 


(a) A single card, format (I5,F10.0), containing 
the number of totai load cases (NTLC) that are in the SAP IV 
Output displacement data deck, and a scale factor (SCALEF) 
that is used in scaling the displacement data, must be input 
berore tne displacement deck. This card controls the input 
of SAP IV displacement data tnrough subroutine DATA3. For 
example, a displacement deck from an execution of SAP IV 
that contains the maximum of five load cases would have NTLC 
= 5. The scale factor could be any value desired by the 
user, with default equal to 1. The DATA9 control card for 
tne above case would be a5 in card column five and the 


desired scale factor (SCALEF) in card columns 6-15. 


(b) The deck of dispiacement data obtained from SAP 
IV follows the DATA 9 control card. Since the maxianun 
Number or load cases that are punched out by SAP IV is five, 
the maximum NTLC is five and the actuai number aust be 
specified. Note that if NUDISP, NVDISP, and NWDISP are all 
specified as zeroes through default or in Namelist OPTION, 
the displacement data deck and the DATA 9 control card are 
not required. This feature enables preprocessing only of a 


given finite element model by PSAP. 


(c) The parameter DMAGS in Namelist PICT also 
provides for magnification of displacements. (See section 
ETT .B 9} 

KDATA = 1 
Cails subroutine DArAl, which is prepared by the user 
to read displacement data from a program other than SAP IV 
to be plotted. 
KDATA = 5 
Calls subroutine DATA 5, which is prepared by the user 
to read displacement data from a program other than SAP IV 
to be plotted. 


Use of KDATA =1 or 5 requires that SUBROUTINE PSAP 
be modified to fit the format that the user's input 
displacenent data specifically reguires. A method - for 
accomplishing this change will be discussed in parapgraph Cc. 


of this section. 


The format for Namelist PICT is the same as that 
required for Nanmelist OPTION. It requires a Single card, 
6PICT, followed by the specified parameter cards and the 
&END cards, all cards beginning in card column 2. This 
namelist contains the values needed to specify the type of 
plot that is desired and what information is to be included 
On the flots. A detailed summary of Namelist PIcT is 
contained in Reference 2 and is given here for user 
convenience. 


FORTRAN name - Default value - Description 
KHORZ =~ 1 Integer designating the i 
1 


(see Figure 5.) 


ANWId ONIMAIA Z - X NI NMOHS ‘*COUd IIHAVADONLYO 
ANOITIGO N¥ NOI SNOILVLON ANW WALSAS ALYNIGHOOD - G arnbty 


KVERT = 2 


PHT - 0.0 


THETA - 0.0 


PSI = 0.0 
NEWFR - 1 
ISCALE - 1 


PLOTSZ - 10.0 


XORGN - 0.0 


YORGN - 0.0 


PSCALE - 1.0 


NOTAT - 0 


ALET = 0.15 


KDiSP = 0 


Integer designating the vertical axes of 
the viewing plane where 1=X,, 2=¥,,and 3=2, . 
(see Figure 5.) 
Angular rotation of model about its X-axis in 
degrees (must be performed third, see Figure 
5.) 
Angular rotation of nodel about its Y-axis in 
degrees (must be performed second, see Figure 
3.) 
Angular rotation of model about its Z~axis in 
degrees (must be performed first, see Figure 
3+) 

1 Frame change before plotting 

{a frame change resets the Y-origin past 
previous plot by XSPACE given in Namelist 
OPTION and resets the X-origin at 0.0) 
0 no frame change before plotting 

1 automatic scaling of plot and computation 
of proper origin location 
2 user-specified origin and scaling 
Maximum dimension desired on completed plot, 
in inches (used for scaling if ISCALE=1) 
X-location of plot origin 

(used if ISCALE=2) 
Y-Location of plot origin 

(used if ISCALE=2) 

Model size reduction factor 

(i.e., PSCALE is equal to actual aodel size 
divided by desired plot size, used 

if ISCALE=2) 

0 no numbering on plots 

1 numbering of grid points 

2 numbering of elements 

Height of integers specified by NOTAT, in 
inches (must be 2 0.07) 

0 plot of undeformed structure 
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plot of deformed structure 


1 
2 exploded plot 

3 displacements represented by vectors 
1 


IDMAG ~- 2 direct magnification of displacement data 
by DMAGS 
2 scaling of displacement data to a 
maximum value ofDMAGS 

DMAGS ~- 1.0 Magnification of displacements 
(if KDISP= 1 or 3) 
Reduction factor of elements 
(if KDIS P=2) 

KSYMXY - 0 1 symmetry about X-Y plane 

KSYMXZ - 0 1 symmetry about X-Z plane 

KSYMYZ - 0 1 symmetry about Y-Z plane 


Syametries are perrormed consecutively (i.e., a plate 
quadrant with KSYMXZ and KSYMYZ equal to one would yield a 
complete flate). 


XXMAX,XXMIN - 1.0E+20 Locate cutting planes parallel to 
principal planes 

YYMAX,YYMIN - 1.0E-20 (X-Y, X-Z, Y-Z) to limit plot 

ZZMAX,ZZMIN 


NDMAX - 9999999 Maximum grid point identification 
Humber to be included in the plot 

NDMIN - O Minimum grid point identification 
number to be included in the plot 

NELMAZ - 9999999 Maximum element identification 
number to be included in the plot 

NELMIN - QO Minimum element identification 
number to be included in the plot 

KODE - 0 Specifies control option after a 
plot is complete 


0 last plot, exit from progran 
1 read another Namelist PICT 


2 cead a new set of displacement 
data (see NOTE 1 below) 
3 read a complete new set of input 
data starting with a title card. 
(NOTE 1: For SAP IV displacement data, KODE=2 signifies 
that the next load case displacements will be assigned to 


the model nodal points.) 


The previous sections describe a complete basic set 
of input data, if KOD2 =0 in Namelist PICT. For KODE =1, 2, 
or 3, additional sections of the deck must be repeated. The 
deck must end with a Namelist PICT having a value of KODE = 
Can Ete An example input data deck and output flots for 
the simple truss problem of Appendix A is found in Appendix 
B. 


Cc. METHOD FOR ALTERING SUBROUTINE PSAP 


In the event a user has geometry and displacement data 
decks from a program other than SAP IV, it is possible to 
plot those decks with PSAP. The subroutine will handle 
rod-like elements, triangular elements, or quadrilateral 
elements when they are input in acceptable format. By 
studying subroutines GEON9 and DATA9, and SUBROUTINE PSAP 
(Appenagix D), the necessary sequence of input can be 
determined. PSAP is presently constructed so that a user 
may supply his own routines through the use of subroutines 
GEOM1 or GEOM2 for geometry data and DATA1 or DATAS for 
displacement data. In order to add a subroutine to PSAP, it 
need only be placed after the main calling program in the 
sequence of control cards as discussed in paragraph II.8.2. 
The essential features for the input of the geometry data 
are the nodal points, with their X,Y,Z coordinates, and the 


connectitity sequence for the finite element model. The 


necessary part of displacement data input is the node point 


number, with the U,V,W displacements. Adding a user-prepared 
subroutine through GEOM1 or 2 and DATA1 or 5 should prove to 
be a relatively straight-forward task for the user who 
desires to do so. The listing of GEOM9 in Appendix D can be 


used as a guide. 


De SIGNIFICANT ASPECTS OF SUBROUTINE PSAP 


Often the absence or presence of elements ina 
finite element nodel cannot be determined from a 
conventional obligue orthographic projection. For example, a 
line element that is coincident with an edge of a triangular 
or quadrilateral element could not be detected. Te show 
clearly each element, PSAP contains an algorithna for 
generating exploded oblique orthographic projections. This 


can be a valuable tool in checking the toplogy of an 
analytical model. (i.e., KDISP = 2) 


2. Portions of Models 


The ability to isolate a portion of a model for detailed 
examination is a very useful and desirable asset of a 
preprocessor. SUBROUTINE PSAP has the capability of 
specifying cutting planes (i.e.,XXMIN ,XXMAX ,YYMIN,YYMAX, 
etc.) or maximum and minumum element numbers (i.e., NELMIN, 


NELMAX). Examples of this are shown in Appendix D. 


The specification of view of a model is done through 
the use of the parameters KVERT,KHORZ, PSI, THETA, and PHI 
that are found in Namelist PICT (Section III.8B.9) The 
specific details of how this is accomplished can be found in 
Reference 2. There are a great humber of possible 
combinations of the above parameters, and Appendix C 
illustrates various combinations that were used on a section 
of a finite element wing model. The Calcomp pilots, along 


with the parameters as specified, are found in each figure. 


(a) The most probable source of error could be 
incorrect data deck preparation or deck sequencing. It is 
important to use Figure 2 as a guide while preparing the 


input for PSAP. 


(6) Other errors nay occur if the arrays ZZ and 
DISPD are not dimensioned correctly as discussed in 
paragraph III.B.3. During the execution of the program 
several manipulations are performed with the two arrays, and 
it is possible for addresses of tae data to be lost within 
the IBM 360/67. The 360 System error messages may not 
directly indicate that the dimensions of the arrays have 


been exceeded. 


(c) A most important point to remember in generating 
a sequence of plots is that once a parameter has been 
assigned a value in a namelist, it retains that value until 
it is reassigned. For example, if PLOTSZ is assigned a value 


Of GeO TOE the £lrst. if a series of plots and it is not 


redefined in any subseguent Namelists PICT , the value of 
PLOTSZ will retain the value (8.0), as originally specified. 


(dq) It is possible to make any number of errors, 
however, and all of them cannot be anticipated. The error 
messages from the IBM 360/67 System will, in most cases, be 
straight-forward and facilitate easy trouble-shooting. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


= SS ee ee ee = Se = 


Generali purpose finite element structural analysis 
programs like SAP IV are significant analytical tools that 
can be used most effectively when coupled with a flexible 
pre- and postprocessor. The availability of a powerful 
structural analysis tool with partial pre- and 
postprocessing capability has now been provided for Naval 
Postgraduate School Students. However, at present this 
capability is limited to pre- and postprocessing the data of 
only five of the SAP IV elements. It does not appear 
feasible, at this time, to add any more capability as far as 
the nuaber of SAP IV elements goes. dowever, other 
additional improvements to postprocessing at NPS should be 
made. For example, Appendix C of Reference 2 lists a 
program developed to produce contour plots of stress data 
from finite element models. The program, in all probability, 
could be adapted to SAP IV and the NPS IBM 360/67. This 
represents a possible avenue to further the present 
structural analysis capability for NPS students. SaP IV 
itself can be expanded. There is a great deal of work being 
done in the area of composite materials, and SAP IV 
possesses only very limited orthotropic material capability. 
However, there are subroutines within SAP IV that could be 
modified to reflect more current methods in handling 


ccmposite materials. 


In conclusion, the finite element method in structural 
engineering is the most powerful analysis tool available 
today and needs to be exploited by NPS students in both 
classroom and research work. Programs like SAP IV are 
available at many Navy laboratories, development centers, 
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and aerospace companies throughout the country. The current 


trends in structural 


analysis are toward computer-aided 


technigues, and exposure to a general purpose package for 


students enrolled in 


asset. 


AE 4101 and AE 4102 is a significant 
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APPENDIX A 


EXAMPLE PROBLEM USING SAP IV 


The schematic given below illustrates the example 
problem for which tne data on the following page were 


prepared. 


PIN- JOINTED 
TRUSS 


1 = 100 inches 
E = 10.12+06 psi 
ALPHA = 12. 62-06 psi 


A = 0.5 Sq.sin. 
A 
A 


1.0 sq.in. 
Q./5 SQ.in- 
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APPENDIX 6B 


EXAMPLE PROBLEM USING SUBROUTINE PSAP 


The pin-jointed truss problem in Appendix A is 
as an example in iliustrating SUBROUTINE PSAP 
entire computer card deck used to generate Figures 
8 is listed on the following two pages of this 
The printed computer output for this problem is 


also used 
usé. The 
o, 7, and 
appendix. 


shown in 


Figure 9. Note that the displacements of node number 4 have 


been magnified by 1000. 
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Figure o = DISPLACEMENTS‘NODE4) POR TRUSS MODEL PLOTTED 
AS A VECTOR 
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SAMPLE COMPUTER PRINTOUT FROM PSAP OF TRUSS 
MODEL 


Figure 


APPENDIX C 


SAMPLE SUBROUTINE PSAP OUTPUTS FOR A FINITE ELEMENT MODEL OF 
A WING 


The figures contained in this appendix are of a finite 
element wing model developed in Flight Vehicle Structural 
Analysis II (AE 4102, QTR.IV, 1976). The original wing was 
designed in AE 4274, which is the latter half of the 
Subsonic structural design sequence in the Department of 
Aeronautics. The finite element model of tne wing was 
developed in AE 4102 as a class project. (The structural 
analysis and design sequences in the the Department of 
Aeronautics are constructed so that they interelate the two 
different areas.) Each student was assigned responsibility 
for a portion of tne wing and the student whose design was 


used was the analysis leader. 


Figure 10 shows the axis orientation used in the 
original definition of the wing as designed in Az 4274. 
Figure 11 and 12 show the entire model with over four 
hundred aembrane elements and with different views of the 
aodei specified. In Figures 13 and 14 a single portion of 
the wing is shown with different rotation angles. The 


significant Namelist PICT values used to orient the nodel 


are included in the figures. 


MIG OEP 


WING ROOT 


Figure 10 - WING MODEL AXIS ORIENTATION 
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KHORZ=1, KVERT=3 , PHI=79, 
THETA=-25 , PSI=-10 


Pigure 11 - ORTHOGRAPHIC PROJECTION OF WING MODEL -1 
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APPENDIX D 


SUBROUTINE PSAP LISTING 
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